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Abstract 

The history and phenomenology of hadronic parity nonconserva- 
tion (PNC) is reviewed. We discuss the current status of the experi- 
mental tests and theory. We describe a re-analysis of the asymmetry 
for p+p that, when combined with other experimental constraints and 
with a recent lattice QCD calculation of the weak pion-nucleon cou- 
pling /i^, reveals a much more consistent pattern of PNC couplings. 
In particular, isoscalar coupling strengths are similar to but some- 
what larger than the "best value" estimate of Donoghue, Desplan- 
ques, and Holstein, while both lattice QCD and experiment indicate 
a suppressed h^. We discuss the relationship between meson-exchange 
models of hadronic PNC and formulations based on effective theory, 
stressing their general compatibility as well as the challenge presented 
to theory by experiment, as several of the most precise measurements 
involve significant momentum scales. Future directions are proposed. 



1 Introduction 



The experimental study of parity nonconservation (PNC) in AS* = hadronic 
interactions began only a year after Lee and Yang pointed out that parity 
might be violated in the weak interaction [1]. In 1957, the year that the 
^°Co beta decay experiment of Wu et al. confirmed the Lee and Yang hy- 
pothesis [2], Neil Tanner published the (negative) results of an experiment 
seeking to measure PNC in the reaction ^^F{p,ay^O [3]. Although the Tan- 
ner experiment lacked the sensitivity required to observe a signal, his effort 
was followed by others, with experiments designed to isolate tiny weak ef- 
fects within systems where strong and electromagnetic interactions dominate 
continuing until today. The presence of weak effects is apparent from results 
such as the 2% photon asymmetry in the electromagnetic decay of an isomer 
of isoHf [4] 

A^C^'^m* ^^^m + 7) = -(L66 ± 0.18) X 10-2 (1) 

or the nearly 10% asymmetry in the scattering of longitudinally polarized 
neutrons from ^^^La [5] 

Ah{n + ^3^La) = (9.55 ± 0.35) x 10"^ (2) 

The asymmetries given in Eqs. (1) and (2) are, however, anomalously large, 
amplified by nearly degenerate nuclear states having the same spin but op- 
posite parity. Indeed the natural scale of hadronic PV effects is five orders of 
magnitude smaller, ~ GpF^ ~ 10^'', where F,^ ~ 92.4 MeV is the pion decay 
constant. In the more than half century since the discovery of PNC, many 
searches for hadronic PNC effects have been performed, both with and with- 
out the use of nuclear amplification, creating a substantial body of results. 
(See, for example, the reviews of [6-8]). 

The focus of such work has changed from detection to analyzing the 
structure of the PNC response. At the time of the first experiments, only 
the charged weak current (CC) was known to exist, 

J+ = cos OcU'Jij.il + ^bjd + sin 6'cM7m(1 + 75)s (3) 
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where 6c ~ 13° is the Cabibbo angle and 75 = —i'y^^^j^^^ is the negative of 
that defined by Bjorken and Drell [9]. The AS* = weak interaction resulting 
from the contraction of this current and its Hermitian conjugate, 

?^S^ = ^j;tj+.+h.c. (4) 

then consists of two components, a A7 = 0,2 piece multiphed by cos^ ~ 1 
which results from the symmetric product of isotopic spin 1 = 1 parts of the 
weak current and a AJ = 1 piece multiplied by sin^ 6c ~ 0.04 which results 
from the symmetric product of / = 1/2 parts of the current. Thus the 
A/ = 1 piece of the strangeness conserving PNC weak interaction generated 
by the current of Eq. (3) is strongly suppressed with respect to its A7 = 0, 2 
counterparts. Early on it was recognized that experimental studies of this 
AJ = 1 component could be a probe possible neutral weak currents (NC). 
The development of the standard model and the experimental discovery of 
the hadronic NC 

Jl = u-i^{l + 75)ii - J7m(1 + l^)d - 4 sin^ J^' (5) 

where is the Weinberg angle and J^"* the electromagnetic current, led 
to an important program of semi-leptonic weak interaction studies of this 
current. In the case of hadronic weak interactions, this current generates a 
Hamiltonian 

<^=||^:^^°'^+h.c.. (6) 

that does add an unsuppressed NC contribution to the AS" = 1 interaction: 
the product in Eq. (6) generates A7 = 0,1,2 interactions. As NCs do 
not contribute to flavor-changing interactions, it was recognized that the 
resulting opportunity to study the NC hadronic weak interaction in A^ = 
systems - particularly A7 = 1 channels - is unique. The nucleon-nucleon 
and nuclear systems are the only practical arenas for such studies, with PNC 
then available as a tool for isolating weak effects in systems dominated by 
strong and electromagnetic interactions. Summaries of experimental work 
are given in Sec. 2 [6-8]. 

During the past five decades a great deal of theoretical work has also been 
performed, with the goal of understanding the patterns of observation and 
non-observation of PNC that emerged from the experiments. We describe 
several of the important unresolved issues in our description of PNC in Sec. 



2 



3, using the traditional meson-exchange picture of hadronic weak interac- 
tions. In Sec. 4 discuss issues in the context of partial-wave analyses and 
effective field theory, approaches that allow in principle for a more system- 
atic approach to PNC (though their description of the weak nucleon-nucleon 
(NN) interaction is largely equivalent to the meson-exchange formulation). 
Then in Sec. 5, we describe an envisioned program that would have the po- 
tential to clear up many of the existing uncertainties. Our paper concludes 
with a brief summary in Sec. 6. 



2 Experimental Summary 

In the more than five decades since publication of the Tanner paper, many 
additional experiments have been carried out to probe the hadronic weak 
interaction. Ideally one would do a series of measurements in the NN sys- 
tem, where the strong interaction physics can be handled precisely: in the 
Danilov amplitude decomposition we describe in Sec. 4, a minimum of five 
measurements would be needed to determine the full set of S-P amplitudes. 
As the natural scale of NN PNC observables is ~ 10"''^, exceptional effort 
must be invested to accumulate sufficient counts to see an effect and to con- 
trol systematics that could yield false signals. The results to date include 

a) The analyzing power for the scattering of longitudinally polarized pro- 
tons from an unpolarized proton target has been measured at 13.6 [10] 
and 15 [11] MeV by groups from Bonn and Los Alamos, and at 45 
MeV [12] by a group from PSI. Furthermore, a medium energy mea- 
surement oi Al was made at TRIUMF [13]. The results are 



b) The asymmetry in the radiative capture of cold polarized neutrons 
on a parahydrogen target, np — >■ dj, was measured at Grenoble in 
1977 [14] and much more recently at LANL [15]. The two results are 



^l(pp; 13.6 MeV) 
Al(pp; 15 MeV) 
Ai(pp;45 MeV) 
ylL(pp;221 MeV 



(-0.93 ±0.20 ±0.05) X 10 



-7 



(-1.7 ±0.8) X 10"^ 
(-1.57 ±0.23) X 10"^ 
(+0.84 ±0.34) X 10"^ 



(7) 
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of comparable sensitivity, with both faihng to find a non-zero signal for 
PNC, 

>l7(np)|Grenoble = (0.6 ± 2.1) X 10"^ 

A(np)|LANL = (-1.2 ± 1.9 ± 0.2) X 10-^. (8) 

c) The circular polarization of the 2.22 McV photon emitted in the cap- 
ture of unpolarized thermal neutrons by protons was measured at the 
Leningrad reactor, yielding the upper bound [16, 17] 

P^(np) = (1.8 ± 1.8) X 10"^ (9) 

Because the only definitive measurement of PNC in the NN system comes 
from the pp system, experimentalists have turned to few-body systems. Here 
calculations can be done in principle with quasi-exact nonperturbative strong 
interaction methods, though applications to scattering states are far less de- 
veloped than to bound states. Faddeev-Yakubovsky methods and variational 
methods in hyperspherical harmonic bases can and have been applied to A— A 
systems: the first bench-mark PNC calculations for ■^He(n,p)^H have been 
recently reported [18]. Hyperspherical harmonic and quantum Monte Carlo 
methods are under active development for A=5 PNC applications such as the 
spin rotation of polarized neutrons in helium [18, 19]. Thus within the next 
few years, quasi-exact methods should exist for extracting weak couplings 
from few-body systems. In the interim, model-based scattering calculations 
can be done with bound-state wave functions that are effectively exact. Few- 
body experiments include: 

d) A PSl experiment measured the cross section difference in the scatter- 
ing of longitudinally polarized protons from a "^He target at 46 MeV, 
yielding an asymmetry [20] 

ALipa; 46 MeV) = (-3.3 ± 0.9) x 10"^. (10) 

e) A NIST study of the spin rotation of transversely polarized neutrons 
in passage through a "^He target has yielded the upper bound [21] 

^— = (1.7±9.1± 1.4) X 10"W/m. (11) 
dz 

(As this observable is the isospin mirror of p+^He, it would allow an 
isoscalar/isovector separation to be done for A=5.) 
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f ) An upper bound also has been established [22] on the analyzing power 
for longitudinally polarized protons scattering on deuterium 

^^(pd; 15 MeV) = (-0.35 ± 0.85) x 10~^. (12) 

As so few of the attempted experiments in NN and few-body experiments 
succeeded in isolating nonzero effects at the expected ~ 10~^ level, experi- 
menters turned in the late 1970s and 1980s to more complicated "two-level" 
nuclei where chance amphfications of PNC observables could be exploited. 
The origin of this amplification is level mixing: in the presence of a parity 
violating interaction described by the Hamiltonian "H^ = "H^*" + "H^*", the 
J'^ strong eigenstates of same angular momentum and but opposite parity 
will mix via the weak interaction. The PNC admixture can be determined 
from first-order perturbation theory 

iV'J+i > - 
iV'J-i > - 

where i and k represent all other quantum labels of the nuclear levels. For 
typical nuclear weak matrix elements < 0j-fc|'Hw|0j+j >~ 1 eV, which can 
be compared to the spacings of neighboring opposite-parity major shells ~ 10 
MeV. Consequently the PNC amplitudes in Eq. (13) are typically~ 10~^. 
However, if two states fortuitously form a nearly degenerate parity doublet, 
the mixing can be significantly enhanced. In this case, Eq. (13) also simplifies 
because the two-level mixing will dominate all other contributions. Many 
early experiments utilized nuclei with such parity doublets to compensate 
for experimental sensitivities that would otherwise have been insufficient to 
probe l-Lyj. 

Several experiments yielded nonzero results, and from such two-level sys- 
tems one can often extract the magnitude of the nuclear matrix element of 
T-Lw But the next step - deducing from the many-body matrix element con- 
straints on the underlying weak couplings of T-Lw - can be problematic due to 
the uncertainties in model-based calculations of nuclear wave functions. Yet 
there are possible exceptions: 

g) Four independent efforts were mounted to search for the circular po- 
larization of photons emitted in the decay from the J^, / = 0~, 1.081 



J, Ej+i — Ej-k 



m-i>+Z^ — ^ ^ I0j+fe> (13) 
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MeV excited state of ^^F to the J^,I = 0+,0 ground state. While 
rather stringent hmits were obtained, the groups failed to find a sig- 
nal. The initial state mixes with the nearby ^ / = O"*", 1 level at 1.042 
MeV, so that = 39 KeV. Thus this system isolates the A/ = 1 PNC 
mixing of greatest interest. The results from the four experiments are 
in good agreement [23-26] 

(-7 ± 20) X 10"^ CalTech/Seattle 

. (3±6)xl0-^ Florence , . 

(-10 ± 18) X 10-^ Mainz ^ ' 

(2 ± 6) X 10-^ Queens 

The experiments reached past the DDH best-value sensitivity while 
failing to detect a nonzero result. Nuclear physics uncertainties in this 
system are unusually modest due to a relationship [27] between the 
mixing matrix element for the long-range pion-exchange contribution 
to PNC and the measured axial-charge /5-decay rate for the transition 
from the J^, I Mi = 0+, 11 ground-state of ^^Ne, the isobaric analog of 
the 1.042 MeV state in ^^F, to the 1.081 MeV state in ^^F. 

h) The photon asymmetry in the radiative decay of the polarized J^,I = 
I , I 110 keV first excited state of -"^^F to the J^, I — ^ , | ground state 
was measured in two independent experiments. The mixing here is be- 
tween these two states, so AE =110 keV. The measured asymmetries 
are [28, 29] 

4 _/ (-8.5 ± 2.6) X 10-5 Seattle , . 

^ ~ \ (-6.8 ± 1.8) X 10-5 Mainz ^ > 

This mixing matrix element is a sum of isoscalar and isovector ampli- 
tudes. The isovector contribution can again be related to the measured 
analog axial-charge (5 decay rate of the ^^Ne ground state to the 110 
keV excited state in ^^F. The isoscalar contribution is assumed to scale 
similarly. While there are arguments to support this procedure, the 
axial-charge beta decay argument is on less firm ground for ^^F than 
for i«F. 



i) Two independent experiments looked for circular polarization of pho- 
tons emitted from the radiative decay of the J^, / = | , | 2.789 MeV 
excited state of ^^Ne to the J^,I = |^, ^ ground state. The mixing 
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here is with the nearby 2.795 MeV J^,I — ^ , | excited state with 
AE — 5.7 KeV. The measured circular polarizations are [30, 31] 




(24 ± 24) X 10-^ 
(3 ± 16) X 10-^ 



Seattle / ChalkRiver 
ChalkRiver / Seattle 



(16) 



While this nucleus is sufficiently light that rather sophisticated shell- 
model calculations can be performed, there are indirect arguments 
based on discrepancies between these calculations and measured El 
transitions to/from the states of interest that nuclear structure uncer- 
tainties are not under control [6]. For this reason results from ^^Ne are 
frequently omitted from global analyses of hadronic PNC. 

Finally, a class of novel experiments in heavy atoms has been included 
in some recent analyses of PNC experiments. The nuclear anapole moment, 
a weak radiative correction, generically becomes the dominant V(electron)- 
A (nuclear) coupling in heavy atoms due to its growth with the mass number 
(cx A^^^) and to the relative weakness of the competing tree-level contribution 
from direct Zq exchange [32,33]. The resulting electron-nucleus interaction 
has the form 



where J and p(r) denote the nuclear spin and density and (Jg is the usual 
Dirac operator for atomic electrons, and where the anapole moment domi- 
nates Ktot in heavy systems, because of the A^^^ growth in its contribution 
to Kfot- Nuclear anapole moments can be isolated through the dependence 
of the atomic PNC signal on nuclear spin (e.g., by examining variations in 
PNC signals for different hyperfine levels). As the largest contribution to 
nuclear anapole moments comes from PNC ground-state parity admixtures 
generated by Hw, in principle they provide another test of the hadronic weak 
interaction. Anapole results include 



Hw = —^i^tot cTe • J p{r) 



(17) 



j) limits on the moment of "^^^TX [34,35], 




Seattle 
Oxford 



(18) 



k) and a measurement of the moment of ^^^Cs [36] , 

/«tot(^^^Cs) = 0.112 ± 0.016 Boulder. 



(19) 
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A rather complex nuclear polarizability arises in the theoretical treatment 
of anapolc moments [37-39]. As in the case of the heavy-nucleus results of 
Eqs. (1) and (2), the calculations are based on nuclear models that can only 
partially capture the relevant physics. The accuracy of the constraints on the 
NN PNC interaction derived from the experiments is consequently somewhat 
difficult to assess. 



3 Meson-Exchange Approach 

For the past thirty years the most common theoretical approach in analyz- 
ing PNC experiments has been that based on the meson-exchange model of 
Desplanques, Donoghue, and Holstein (DDH), developed in 1980 [40]. DDH 
constructed a PNC potential based on tt^, p, and ui exchanges with strong 
vertices described by the Hamiltonian 

= iQ.nnN^^t -TiN + QpN (^^^ + i^a^,k^')f ■ p'^N (20) 
+ 9.N(^^^ + i^a,,k^)u^N (21) 

They chose for their strong couplings gl^j^f/ifr ~ 14.4 and gp/4:7i = ^g'^/iir ~ 
0.62 and, using vector dominance to connect with the electromagnetic inter- 
action, Xv — K'p — K'n — 3.70 and Xs — K'p + K'n — —0.12. 

The weak vertices are based on a Hamiltonian with seven phenomenolog- 
ical couplings 

+ N (^h^T ■ r + + ^ (3r,p^ - f • pn^ 7m75A^ 
+ N {hluj>^ + hlr^uj^) 7^75A^ - h};N{T x p^^)f-^^,N. 

(22) 

Couplings to the neutral pseudoscalar mesons tt", rf , rf' are absent due to 
Barton's theorem, which requires such couplings to be CP violating [41]. 
Combining Eqs. (21) and (22) and doing a Fourier transform yields the 
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coordinate-space DDH potential, 



V2 



(ai + (T2) 



9p ( h% • ^2 + /i! 



1 / + r2 



+ i(l + Xy)ai X (T2 



2 

2m N 

Pi - V2 

1 / n + r2 ^ 



Pi -P2 
(3Tf T| - fi • T2) 



2V6 



,Wp{r) 



X \^(ai - <T2) 
+ i{l + xs)^i x ^2 ■ 



Pi -P2 
2m M 
Pi -P2 
2mN 



+ 



Tl - T2 



-9php i 



ll. riXT2 



{^1 + ^2) ■ 



2m N 

Pi -P2 
2m N 



,Wp{r) 



where Wi{r) — exp(— mir)/47rr is the usual Yukawa potential, r — \fi — 
is the relative NN coordinate, and Pi — —iVi. The resulting PNC NN 
interaction is described in terms of seven phenomenological weak NN-meson 
couplings — /i^, hy, h^^' — though the constant h^' is generally not included, 
since it is a short ranged piece of the dominant pion coupling and a simple bag 
model estimate has shown that it is small [42] . This potential is very closely 
related to model-independent, threshold S — P interactions, as discussed 
below. 

The results of PNC analyses are often presented as in Fig. 1, in terms 
of constraints on the weak couplings. However, the weak NN amplitudes 
depend on a product of weak and strong couplings, so that comparisons of 
extracted weak couplings can become problematic if experiments are not 
analyzed with a common set of strong coefficients. We will return to this 
point in later discussions. 
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DDH [40] 


DDH [40] 


DZ [43] 


FCDH [44] 


Coupling 


Reasonable Range 


"Best" Value 






K 


^ 30 


+ 12 


+3 


+7 


K 


30 ^ -81 


-30 


-22 


-10 


K 


-1 


-0.5 


+1 


-1 


K 


-20 ^ -29 


-25 


-18 


-18 


hi 


15 ^ -27 


-5 


-10 


-13 


hi 


-5 ^ -2 


-3 


-6 


-6 



Table 1: Weak NN-meson couplings as calculated in Refs. [40,43,44], in units 
of GFFll1\f2 ~ 0.38 X 10"'^, where Gp is Fermi's weak coupling and F^^ the 
pion decay constant. 



Using quark model and symmetry methods, DDH attempted to calculate 
values for these weak parameters. However, because of strong interaction un- 
certainties, an accurate estimate proved impossible, so DDH instead quoted 
reasonable ranges and "best values" - or perhaps more accurately, "best 
guesses" - for each. These are listed in Table 1 together with values esti- 
mated by some other investigators. One issue relevant to later discussions 
is that the DDH best value for h}^ is higher than found by some others. 
Skyrme and quark-soliton models tend to give smaller values for /i^ [45,46]. 
A large- Ai'c result has also been obtained recently [47]. 

Analysis of key experiments described in the previous section in terms 
of these unknown weak couplings has produced the well known graph [37] 
shown in Fig. 1, which exploits the fact that certain isoscalar and isovector 
combinations of the DDH parameters dominate the experimental observables. 
This plot dates from 2001 and as we will show later, should be modified to 
reflect more recent developments in PNC studies. 
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-5 

-2 2 4 6 8 10 12 14 

h^^- 0.12 -0.18 

Figure 1: Experimental constraints on linear combinations of isoscalar and 
isovector DDH couplings (in units of 10~^), taken from the 2001 work of [37], 
displaying bounds from four experiments where it is believed that theoretical 
analysis uncertainties are under reasonable control: pp, pa, ^^F, and ^^F. The 
small shaded triangle is consistent with all four experiments. The DDH best 
value point is also shown. Later we show that the data on p + p subsequently 
obtained at TRIUMF [13] and the analysis of Ref. [54] have a significant 
impact on this plot. 
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For the reasons described previously, the most rehable constraints on 
PNC come from a hmited number of measurements where PNC has been 
seen or sharply limited, and where the analysis of the observables can be 
done reliably either with ab initio techniques or by relating observables to 
other measured quantities such as axial-charge beta decay rates. Figure 1 - 
a figure that has been used frequently over the past decade, despite its 2001 
vintage [37] - shows that the allowed bands for these observables - pp, pa, 
^^F, and ^^F - intersect in a limited region defined by linear combinations of 
isoscalar and isovector weak couplings. The resulting value of the isoscalar 
coupling + 0.7hPp) ~ (15 - 18) x 10~^ can be compared to the DDH 
best value ~ 12.5 x 10"''' and reasonable range of (-15.4 -H- 39.0) xlO'''. The 
two alternative sets of theoretical weak couplings given in Table 3 yield ~ 
11 x IQ-^ [43] and ~ 7.3 x IQ-^ [44]. 

If we consider the extended region where at least three of the four bands 
overlap, this range would expand to ~ (13 — 25) x 10"^. Thus the conclusion 
from the 2001 analysis is that the experimental value is likely comparable to 
or somewhat larger than the DDH best value, though quite consistent with 
the reasonable range. One also sees that there is some tension between the 
isoscalar coupling derived from p + p asymmetry and that derived from other 
observables: the p + p asymmetry tends to favor somewhat weaker isoscalar 
couplings. This tension will be re-examined in the next section, in light of 
results on p -|- p obtained after the analysis of [37] that led to Fig. 1. In 
deriving the p + p band, the experimental result was "marginalized" over 
degrees of freedom not shown in the plot, resulting in a broadening of the 
band: this includes the isotensor contribution to p + p as well as a small 
isoscalar contribution orthogonal to the combination — (/i^ -l- 0.7/i^). New 
results discussed in the next section allow improvements in this process. 

The discrepancy in the case of the long-range isovector pion coupling hi. 
is much more substantial. The range of values consistent with the combined 
results of the four ^^F experiments, = (1.2 ± 3.9) x 10"^ [26], and the 
relationship established by measured P decay rates with the inclusion of all 
identified theory uncertainties, \P^\ = (4.2 ± 1.0) x 10^ \hl\ [28], yields 

\hl\ < 1.3 X 10-^. (24) 

Thus \hl\ is substantially smaller than the DDH best value estimate, 4.6 x 
10^^. Although ^^F is a complex nucleus, it is difficult to attribute the dis- 
crepancy to nuclear structure uncertainties because the PNC mixing matrix 
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element needed in the analysis can be determined from the analog beta de- 
cay of ^'^Ne [6,27]: with this constraint, the residual model dependence - 
that is, the variation that can be achieved by using any of the wide range 
of ^^F wave functions to relate ^^F PNC to ^^N (3 decay - was found to be 
< 14% [6]. The axial-charge /3 decay argument relates bare couplings, and 
thus eliminates concerns about operator renormalization that accompanies 
shell-model calculations. 

This hi result is supported by two others. A lattice QCD calculation of 
/i^ was recently performed [48] for a lattice size L ~2.5 fm, lattice spacing 
Qs ~ 0.123 fm, and a pion mass ~ 389 MeV, accomplishing a goal 
envisioned some time ago [49] . The result obtained 

/i^ (connected) = (1.099 ± 0.505 (stat) (sys)) x 10"^ (25) 

is consistent with the ^^F result and also well below the DDH best value. This 
calculation is a first step: it did not include nonperturbative renormalization 
of the bare parity- violating operators, a chiral extrapolation to the physical 
pion mass, or contributions from disconnected (quark loop) diagrams, though 
it was argued that these omitted effects would be expected to lie within the 
designated systematic error. Nevertheless, the result should be considered 
preliminary, pending future work with physical pions and including discon- 
nected diagrams 

The ^®F result is also consistent with the interpretation of a recent mea- 
surement of a PNC triton emission asymmetry coefficient in the reaction 
^Li(n,Q;)^H with polarized cold neutrons. From the nonzero result obtained, 
opNC — (—8.8 ± 2.1) X 10~® [50], the following hmit was obtained 

\hl\ < 1.1 X 10-^. (26) 

However, we are reluctant to adopt this result as independent evidence for a 
suppressed h]^ because the reaction theory on which the limit is based is rather 
simple in its use of a cluster decomposition and schematic interactions. The 
experimental result provides good motivation for further theoretical work. 

4 The TRIUMF 221 MeV p+p Measurement 

A somewhat discouraging aspect of the effort to understand hadronic parity 
PNC has been the slow pace of results since the 1980s. The decade of the 
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1980s gave us several important measurements: the longitudinal asymmetries 
for p+p and p + a, the significant limits from the Queens and Florence 
group on P^(^^F) that pointed to a problem with /i^, our cleanest handle 
on neutral current contributions to PNC, and the measured for ^^F. But 
the last twenty years have proven more difficult. A promising LANL effort 
on Aj{np) - an experiment that could confirm conclusions drawn from ^^F 

- in the end yielded a limit comparable to that established in the Grenoble 
effort of 1977. Another effort is now underway [51] at the SNS cold neutron 
beamline that hopefully will meet with greater success. The NIST experiment 
on neutron spin rotation in He also provided only an upper bound: had 
PNC been detected, a separation of isoscalar and isovector contributions to 
PNC could have been made by combining this result with the existing p+^He 
measurement. As significant difficulties were encountered in this experiment, 
new ideas may be needed before another attempt can be made [52]. The 
anapole moment of Cs was measured - a remarkable feat that does establish 
the expected sharp growth in this weak radiative correction with increasing A 

- but it is difficult to envision any confident determination of PNC couplings 
from the measurement given the complexity of the nuclear polarizability that 
governs this moment. Finally, the recent measurement of PNC in ^Li(n,Q;)'^H 
is a significant advance, but ab initio methods for handling the associated 
analysis are not yet in hand. 

One significant new result is the TRIUMF medium-energy measurement 
of p + p, reported in 2001 and 2003 [13,57], 

Al{221 MeV) = (0.84 ± 0.29(stat) ± 0.17(syst)) x 10"^. (27) 

This measurement was performed at an energy where the ^Sq — '^Pq amplitude 
nearly vanishes, and consequently where the ^P2 — ^D2 partial wave dominates 
the scattering. We discuss this result in detail here because, as explained 
below, this result has been combined with others in a way that has obscured 
its impact on Fig. 1. 

The result, together with the lower energy p + p measurements, was ana- 
lyzed by Carlson et al. [54] using several modern strong potentials - Argonne 
vl8 (AVIS) , Bonn-2000, Nijmegen-I. The full scattering problem with strong. 
Coulomb, and weak potentials was solved in both coordinate and momentum 
space. This work remains state-of-the-art, the definitive treatment of PNC 
in p + p. The approach resembles early work by Driscoll and Miller [53] and 
confirmed many of their conclusions, including the importance of distorted 
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waves even at rather modest momentum transfers [53]. While channels up 

to angular momentum J=8 were included, it was shown that the region up 
to the laboratory energy of 221 MeV could be described accurately by re- 
taining the first two PNC partial waves of the DDH potential, — ^Pq and 
'P2 - 'D,. 

The results of Carlson et al. were presented as bounds on DDH weak cou- 
plings, and those bounds appears roughly consistent with the pre-TRIUMF 
analysis shown in Fig. 1. For this reason Fig. 1 is frequently used today 
as a summary of PNC constraints. However the analysis of Ref. [54] used 
CD Bonn strong couphngs that differ substantially from those of the DDH 
potential, 

„CD Bonn CD Bonn ^CD Bonn ^CD Bonn 

ZP 1 16 ^ 1 89 ^ 1 65 — = (281 

r,DDH " nDDH °^ ^DDH ^'^^ ^DDH V^°; 

9p Qui Xv Xs 

This, of course, does not affect the validity of the calculations - only their 
interpretation when combined with calculations that may have made other 
choices in their strong couplings. The fact that the Carlson et al. obtained 
numerical values for weak couplings that were roughly compatible with those 
shown in Fig. 1 is likely the reason Fig. 1 remains in wide use. In fact, when 
the differences in strong couplings arc taken into account, the TRIUMF data 
and the Carlson et al. analysis require an important revision of Fig. 1 
that brings the experiments into better agreement, as we explain below. (We 
note that the confusion that can occur when comparisons of experiments and 
calculations are done at the level of weak couplings is an old issue. To avoid 
such confusion, the authors of [6] advocated quoting experimental constraints 
in terms of the product of weak and strong couplings, but this suggestion did 
not catch on. As pionless effective field theory treatments essentially force 
comparisons to be made at the NN amplitude level, not at the weak vertex 
level, the situation may be changing, as we discuss in Sec. 6.) 

The overall strengths of the p and u couplings that enter into p + p 
scattering depend on the isospin combinations 

9,h^ = 9,{hl + + ^) ^'J = aM + hi) (29) 

These parameters appear in the DDH potential alone as well as in the com- 
binations 

QphfXv gXxs. (30) 
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Apd(x|07) haoPP(x|00 

Figure 2: Left panel: 68% and 90% c.l. constraints on the ^5*0 — ^Pq and 
^Pa - coefficients Asp = hpgp{2 + Xv) + K9w{2 + xs) and ApD = 
hpQpXv + h^QujXs derived from the 13.6, 45, and 221 MeV p+p measurements 
and the calculations of [54]. Right panel: Constraints on the DDH potential 
weak couplings imposed by the results on the left when interpreted in terms 
of DDH strong couplings. 
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Thus the parameter space has four degrees of freedom. However in the plane- 
wave Born approximation this collapses to two degrees of freedom associated 
with the coupling combinations 

Asp = gphf{2 + xv) + 9.fC{2 + Xs) Apd = Qphfxv + g.fCxs, (31) 

the coefficients of the — ^Pq and — ^£^2 amplitudes. Thus in this 
limit it is possible to define strictly equivalent sets of weak couplings for 
arbitrary variations in the strong couplings gp-iXViQuj-iXs- This appears to 
be the procedure Carlson et al. used in defining their DDH-equivalent and 
adjusted-DDH couplings. It is also the approximation we make here. 

We use the ^^'o - ^Pq (J=0) and ^^'o - - ^^2 (J=0+2) resuhs 

of [54] which are given separately in the paper. (The J=0-|-2 result is virtu- 
ally identical to the full result that includes all angular momentum channels 
though J=8, showing that only two channels are important.) Treating the 
experimental data as described in [54] and fitting the three experimental data 
points, we find the constraints on Asp and Ap^ displayed in Fig. 2. This 
procedure thus gives us constraints on the product of weak and strong cou- 
plings, the quantities directly determined by experiment. The procedure was 
cross-checked by plugging in CD-Bonn strong couplings and generating the 
analog of Fig. 8 of [54] (the x^ ellipses for and when CD-Bonn strong 
couplings are adopted). The agreement is quite good. We can then derive 
the similar constraints on and for DDH strong couplings, shown in 
the right panel of Fig. 2. The 68% and 90% confidence level contours are 
based on the appropriate x^ for three data points and two parameters (one 
degree of freedom). 

The semi-major and semi-minor axes of the ellipse provide the orthogonal 
constraints 

0.710/1^^^^-^ + 0.705/1^?'^^^ = -18.63 ±1.90 



0.705/i^f - 0.710/i^f = -40.12 ± 19.55 (32) 

at 68% c.l. While the first constraint is close to that needed for a revised Fig. 
1, wc can combine both results to limit the quantity of interest for a revised 
Fig. 1. That is, because of the 221 MeV TRIUMF result, we can bound any 
needed linear combination of /i^^ and h^. From the axis rotation illustrated 
in Fig. 3 we find 

/if + Q.lh^J = -30.75 ± 4.75 (33) 
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Allowing the small isovector and larger isotensor contributions to this quan- 
tity to vary arbitrarily throughout the DDH reasonable ranges, one finds 

-{hi + 0.7/1^ ^^^) = 25.91^?, (34) 

a result that can then be placed on the revised limits graph shown in Fig. 3. 
The new analysis removes the tension we observed in Fig. 1. The band 

of values for —{h^ -D-D-H" _|_ Q,Th^^^^) consistent with the constraints from 
]} + a, ^^F, and ""^^F, ~ (15 — 31) x 10^ is in excellent agreement with that 
obtained from p in Eq. (34). The region allowed by all four experiments, 
~ (20 — 31) X 10^, is centered on a value that is approximately twice the 
DDH best value. 

Figure 3 also includes Wasem's recent lattice QCD calculation [48] of h]^, 
shown as the blue vertical band. We have plotted this result by adding the 
systematic error to the statistical error, and note that it is difficult to estimate 
the effects of the omitted disconnected diagrams and of the nonphysical pion 
mass employed (389 MeV). 

Even without new experiments, there are important steps that can be 
taken to further tighten the bounds shown in Fig. 3. First, Wasem's bench- 
mark calculation has demonstrated that lattice QCD has reached the point 
that it can impact hadronic PNC analyses. The completion of the hi work 
- the evaluation of the disconnected contribution and the repetition of the 
connected calculation with physical or near-physical pion masses - appears 
to be feasible with computing resources now coming available (though the re- 
quirements may be 10-100 times the CPU investment made in the connected 
calculation, where the contractions required about 6 months of running on 
LLNL's Edge GPU cluster). There is also a second important opportunity: 
the p + p band shown in Fig. 3 requires one to remove the isovector and 
isotensor contributions from /i^^ + 0.7hf^. This was done by allowing these 
contributions to vary over their full DDH reasonable ranges, resulting in a 
substantial expansion of the uncertainty on /i° -|- 0.7/i°. As the isovector 
couphngs are thought to be very small, this step is controlled by the DDH 
reasonable-range uncertainty on /i^. This parameter has no disconnected 
contribution in lattice QCD and thus should be calculable to high accuracy. 
The completion of a lattice QCD calculation of /i^ could thus significantly 
narrow the p + p band in Fig. 3. The CalLat Collaboration has proposed a 
program of lattice QCD calculations focused on hi and h^ [58] . 

Second, we would advocated that the authors of [54] revisit that impor- 
tant work. It should be possible to express the p + p asymmetries calculated 
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Figure 3: Left panel: The 68% and 90% c.l. constraints on h'^^ and 
from the right panel of Fig. 2 are rotated into a form that can be use in 
redrawing Fig. 1. Right panel: Our suggested replacement for Fig 1 using 
the new analysis of p+p to constraint the isoscalar combination /i°+0.7/i° and 
showing a recent lattice QCD estimate of the connected-diagram contribution 
to h\ (blue vertical band). Units are 10"''. The gray shaded area is consistent 
with all experiments. 
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there analytically, as coefficients associated with combinations of weak and 
strong parameters, hpg^. hpQpXv, hi^g^^, and h^g^^xs-, evaluated for each of 
the three energies of interest (13.6, 45, and 221 MeV). The coefficients could 
be tabulated for the three strong potentials explored in [54] and for a range 
of reasonable single-nucleon form-factor masses. (A form factor was used 
in [54] because of the large relative momentum transfer at E'lab = 221 MeV.) 
This would allow one to simplify the analysis described above - and guaran- 
tee its correspondence with the exact numerical results. The dependence of 
the expansion coefficients on the choice of strong potential and on the value 
of the form-factor mass would provide a good measure of strong-interaction 
uncertainties. 



5 The DDH Potential and its S-P Reduction 

At sufficiently low energies in elementary two-nucleon systems the matrix 
elements of the PNC interaction can be written in terms of five elementary SP 
amplitudes, as Danilov emphasized in early work. Because long-range pion 
exchange contributes to the ^5"! — ^ Pi channel, the short-range structure can 
be more readily distinguished in this amplitude, through the contribution of 
higher partial waves. Thus the zero-range approximation to the low-energy 
PNC interaction is sometimes elaborated through the addition of a sixth 
parameter, the pion range. 

As recent treatments based on lowest-order effective field theories (EFTs) 
begin with the threshold behavior of PNC amplitudes, here we take a similar 
approach with the DDH potential, to illustrate the compatibility of the po- 
tential and EFT approaches. By doing a Taylor expansion on the initial and 
final states, the DDH potential can be expanded in a power series in deriva- 
tives compared to m, where m represents m,r, i^p-i or m^;, while the DDH 
Yukawa potential effectively contracts to a contact form. The Danilov limit, 
retention of the lowest order contribution corresponding to S4->P amplitudes, 
corresponds to the first term in the expansion where the identification 

p-rnr -i 



is made in Eq. (23). [Whether one considers the delta function written above 
as a true delta function or the representation of a "fuzzy" one, m^e~"*^/47rr, 
is really a matter of taste [55] , in that PNC will always be treated perturba- 
tively] 
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A bit of work is needed in expanding the DDH potential to lowest order 
because ten operators arise in the reduction of Eq. (23), but only five of 
these are independent in the S-P limit. Defining 

Vs S{f) =V S{f) + 6{f) V Va S{f) = -V S{r) + 6{f) V 

one finds the following identities (valid in S-P matrix elements) 

•f-)- ^ 

V A S{f) ■ {ai - a2)fi ■ vz = - Va 5(f) ■ (<?! - (T2) - 2i(CTi x ^2)- Vs 5(^) 

V s S{f) ■ i{ai X a2)ri ■ r2 = -2 V a 5{f) ■ {Bi - Bz) - i{ai x az)- V s 8{r) 

VA5{r)-{a,-a2){Tl + T^) = -VsS{f^-i{a,xa2){r^ + r^) 

Va 5(f)- (0^1 + ^2) (rf-rl) = -VsS{r)-{a, + a2)i{nXT2), 

Va 5(f) ■ ((Ti - a2){fi (g) f2)2o ^ - Vs 5(f) ■ i{ai x 0^2) (fi ® 7^2)20 (35) 

where the rank-two tensor product (fi (g) 7^2)20 = (3Tfr| — fi • T2)/VQ. 

These identities can be used to write the most general lowest-order (LO) 
PNC potential in many equivalent ways. We choose to retain the five oper- 
ators used in [6], using Eqs. (35) to eliminate the remaining five, obtaining 



35 ip / 1 Va 5(f) I Vs 5(f) ^ 

, A=^Si-^Pi 1 Va 5(f) .^ ^w . 



i5o_3p 1 Va 5(f) , 

A2 °-- 5- • (CTI - (72)(Ti (8) T2)20 



i 2mN rn?p 



(36) 



where the subscripts on the coefficients A denote the change in isospin A/ 
induced by the operator, and the superscripts denote the transitions induced 
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(37) 



by the operators. (The hnear combinations of the two isoscalar operators 
formed in the first two hnes above ensure simple projections.) This hermi- 
tian S-P PNC potential is the most general contact interaction that can be 
constructed, once the identities of Eq. (35) are used to remove redundant 
operators in favor of the five we have chosen. Consequently the 1/m expan- 
sion of the DDH potential that keeps only the LO terms must have this form. 
This leads to an identification of DDH weak couplings with the coefficients 
of the potential, 

^Win' = -5p/^p(2 + Xv)-5./i°(2 + X5) 
K%DH° = -g,hli'2 + xv)-9.hl{2 + Xs) 

K'h-DH" = -9phl{2 + xv) 



On the right we have indicated the DDH "best value" equivalents as a very 
rough guide to magnitudes, with the main points being the strength of the 
long-range ^i]jdh^ coupling (keep in mind that this often contributes only to 
the exchange channel, which can weaken the contribution by a factor ~ 5 [6]) 
and the weakness of ^i^j^^h'- Here we have taken ~ TUp, to streamline 
the expressions: if one does not make this simplification, all cu terms are 
multiplied by rn?p/mlj ~ 0.97. 

Sometimes the DDH potential is discussed as though it were quite distinct 
from the nonrelativistic LO potential of Eq. (36), which can be constructed 
from symmetries and power-counting in external momenta. But the expres- 
sions above show that this is not the case. First, in the isoscalar and isotensor 
channels, there is a 1-to-l mapping of DDH weak couplings and LO operator 
coefficients. 

In the isovector channel in principle there is an overcompleteless in the LO 
reduction of the DDH theory because four weak couplings map onto the two 
independent operators 



DDH 1 ^^l DDH i ^ \"'-Ki "'pj "'UJ1 "'p J ■ 
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However JS^i]^^^^ is completely dominated by h\, with the contributions 
from /i^, and /i^' contributing at the level of a few percent. (In fact 
customarily /i^' = 0.) Consequently, if the DDH potential were used in 
processes where the momentum scale g — )■ 0, one would find that A/ = 1 
observables would depend only on two parameters, h\ and the combination 
h^p + ^i(5'a;(2 + Xs))/{gp{'^ + Xv))' with this observation the DDH potential 
would be functionally equivalent to the most general LO effective theory. 
(The momentum scale for most PNC observables is not insignificant, as we 
will discuss later, so the DDH potential does take partial account oi P — D 
and other contributions, as discussed in the next section.) 

This discussion also further clarifies the motivation for our reanalysis of 
the Carlson et al. calculation of p + p: the mapping from effective operator 
coefficients to DDH parameters involves products of weak and strong coeffi- 
cients. Any consistent global analysis in terms of weak couplings must use a 
fixed set of strong coefficients, so the these products are properly determined. 

The effective Hamiltonian of Eq.(36) is written in terms of operators of 
definite isospin, but one can of course use it to calculate NN observables. For 
initial and final states with strong-interaction distorted waves one finds 



^ mj, I zniN 
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i 2mN 
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where we have defined 



A i5o-3Po 

PP ~ "I" -'^-i "I" 

A i5o-3Po 

A;f°-^^° = Ai^-^^o + A;^°-^^o (39) 

In terms of quantum number \{LS)JMj;TMt), here the pn states have 
been defined as the normahzed states \^So)pn = | (00)00; 10), |^Pi)p„ = 
|(10)lMj;00), etc. 



6 DDH Potential &6 EFTs: A Rosetta Stone 

We have seen, with the customary choice in the DDH potential of h^^ = 0, 
the one redundancy in the potential, when viewed in the S-P limit, is in the 
isovector channel. This is the channel dominated by pion-exchange, and thus 
the channel where the contact-gradient expansion that comes from taking 
the limit mmeson — ^ oo will fail first, as a function of increasing momentum 
transfer. Following [6] , this redundancy can be expressed as an invariance of 
the DDH S-P threshold amplitudes under the simultaneous shifts 

g^NNhl — >■ g-KNNh\ -\- T] 



9phl 9ph\ 



^/2ml 4: + iJ,s + IJ'V 



This transformation leaves the isosector <H-^ Pq and <H-^ Pi amplitudes 
invariant. While the triplet isovector channel contains the pion contribution, 
the S — P limit provides no information on the range of the interaction: the 
pion contribution cannot be disentangled from the p/u contribution. 

The are two P — D AT = 1 partial waves sensitive to the pion-range 
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Here ® denotes a tensor product. These operators come from the next order 
in the Taylor expansion of the initial and final wave functions. The derivation 
of these terms from a potential model would identify the delta functions in 
these terms as the r^- weighted moments of the Yukawa potentials e~"*^/r, 
which would consequently emphasize pion-exchange contributions. 

In principle, these NLO contributions break the degeneracy among the 
isovector parameters of the DDH potential. Given a complete set of low- 
energy data that includes some observable with sensitivity to D-waves but 
where momentum scales still allow a Taylor series approximation, the data 
will effectively fix the strength of gT^MNh\. The "sixth degree of freedom" in 
the DDH potential thus allows one to take account of some average effect 
of the ^Pi — ^Di and — channels: strong interaction effects that 
differentiate these channels, such as spin-orbit interactions, would not be be 
treated in such a fit. This view of the DDH potential has much in common 
so-called hybrid EFT methods, where conventional potential models are used 
to generate wave functions, while interactions are expanded systematically, 
then evaluated between potential-model wave functions [56]. In practice, 
we will argue below that existing data provide perhaps only two significant 
constraints on PNC, and thus fall far short of what is needed for even a LO 
fit, let alone one that tries to address the finite-range effects of ttItt. 

There have been a number of recent treatments of PNC that have devel- 
oped S — P representations of interactions from a "bottom up" approach, 
developing the S — P amplitudes in the framework of pionless EFT, with 
the inclusion of one derivative. The first such effort was by Zhu et. al [59] 
(though we note EFT approaches have connections to earlier work [60,61]). 
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This treatment retained operators that were related under operator identities 

analogous to those given in Eq. (35), and thus is somewhat more difficult to 
use because of the redundancies. The long-wavelength form of the resulting 
PNC Zhu potential is 

T^/o'" = -2^^(n X f2),(ai + a^) ■ ^ V^^lf) + 2^(n • f2){a, - a^) ■ \V Am 

+ 7-2)^(^1 X aa) ■ ^ Vs5(f) + 2\/6^(ri ® r2)2o(^i - ^2) ■ A^if) 
-2v^||(ri ® r2)2o«('?i X as) ■ Vs5(f) + 2^(ai - <?2) ■ t V^^lf) 
-2^^(ai X a,) ■ \Vsm + + - ^2) • ^ V^<5(f) 



+ r2^)^(ai X a,) ■ sSi^) + ^^\,^'\ tI - r^^){a, + 0^2) • J Va5(0 

(41) 



Girlanda [62] then provided a treatment that addressed the necessary 
operator identities, generating a potential with the requisite five low energy 
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constants. The analysis begins with the twelve operators 



O2 - ^7'^75^5"(^r3a^.^) 
Oz = ^ra7''^/'^r"7/,75V^ 

= ^T„7'^75^a^(^TV^,^) 
C4 = '0T37''^'07^75V' 
O4 = V5T37'^75V'5'^(^t^M-V') 

Oe = ieabzi^TaYi^Wii^na^^i^) (42) 
generating the general PV NN Lagrangian 

6 

Y^iGiOi + GiOi]. (43) 
i=i 

Then, with the use of Fierz transformations and the frcc-particle equation of 
motion, six conditions relating these operators were identified 

62-64 = -2mNOe-6e 

6, = O5 (44) 

Finally, using the feature that the operators Oq and 6^ have the same form 
in the lowest order nonrelativisitic expansion, one determines an effective 
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(pionless) Lagrangian, which reduces to the nonrelativistic form 



V'Girlanda 
LO 



-2^il - V5 d{f) ■ i{a, X 0^2) + [2gi] - S{r) ■ {a, - 0^2) 



+ 



20 



(45) 



(As was done by PhiUips, Schindler, and Springer [63] , in Eq. (44) the factor 
of 1/A^ used by Girlanda has been absorbed into the coefficients, making 
them dimensional.) 



Tabic 2: The coefficients of the S-P PNC potential of Eq. (36) in the 
DDH potential, Girlanda, and Zhu descriptions. Note that multiplicative 
factors of 2mjv?«^ and 2mArm^/A^ must be applied to the Girlanda and 
Zhu entries, respectively, to obtain the dimensionless coefficients A, e.g., 
K%DH = 2(e?i + ^i)[2m^m^] = 2(Ci +€1+63 + C3) [2mArm^/A^] . 



Coeff DDH Girlanda Zhu 



^'■O DDH 


-gph^pi2+xv) - 9M2+XS) 


2{Qi+gi) 


2(Ci+Ci+C3+C3) 


^*-0 DDH 


du^hlxs - Sgph^pXv 


2{g,-gi) 


2(Ci-Ci-3C3+3C3) 


\^So-^Po 
^^1 DDH 


-gphl{2+xv) - gM2+xs) 


g2 


(C2+C2+C4+C4) 


■''^1 DDH 




2ge 


(2C6+C2-C4)) 


^^2 DDH 


-9phl{2+xv) 


-2v^6?5 


2v^(C5+C5) 



Returning to the "canonical form" of the S — P contact potential in 
terms of the partial- wave operators of Eq. (36), the relationships between 
the DDH, Girlanda, and Zhu forms of that potential can be summarized in 
terms of coefficients of that potential, as shown in Table 2. In using this 
table it should be remembered that the DDH results include the assumption 
that a one-boson exchange potential operates between strongly interacting 
initial and final nuclear states. There are contributions from crossed-pion 
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diagrams and delta intermediate states that cannot be factored in this way. 
Inclusion of such terms would alter the mapping between coefficients and 
meson couplings shown in the table. See, for example, Refs. [59,64]. 

7 Observables and Momentum Scales 

The decomposition of the PNC interaction into SP amphtudes provides an in- 
teresting way to think about the need for additional experimental constraints 
on the PNC interaction. This type of approach - envisioning experiments 
that might constrain the five degrees of freedom of an S* ^ P potential, valid 
near threshold - dates back to the work on Danilov [61]. (See also Ref. [60].) 

One can approach the problem of obtaining five experimental constraints 
on the PNC potential with varying degrees of reahsm. For example, Danilov 
suggested a treatment where the strong-interaction input would be limited 
to the "^Si and ^5*0 strong phase shifts. The parameters derived from fitting 
experiment in this limit would encode and thus entangle weak and strong 
physics: for example, the strong short-range repulsion that carves out the r ~ 
0.5 fm hole in the nucleon-nucleon correlation function would be entangled 
with the short-range weak physics the DDH potential attributes to p and u 
exchange. Yet the effective couplings derived from fitting data would still 
provide a valid parameterization of low-energy weak NN interactions, and 
once those parameters were determined, could be used to make predictions. 

Several such "unitarized" strong phase-shift methods designed to satisfy 
the generalized Watson theorem were developed and employed in calculations 
of PNC observables such as Al{p + p) in the 1970s and 80s [65,66]. The 
inadequacies of such approaches in applications to low-energy observables 
such as Al{p + p) were noted in the 1980s [6]. The importance of treating 
the strong distortion of the partial waves was first demonstrated in explicit 
calculations by Driscoll and Miller [53]. The need for distorted waves follows 
from simple considerations: our highest precision measurements of Al{p + p) 
were performed at 13.6 and 45 MeV, or at center-of-mass relative momenta 
of 80 and 145 MeV. A treatment like that envisioned by Danilov in which the 
information about the strong interaction is encoded entirely in the scattering 
length cannot be successfully applied to these data: the range of validity is 
determined not by the naive estimate of m^, but by the anomalously large 
scattering lengths of the nn, np, and pp channels of ~ 20 f . Using the effective 
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range expansion, one find that the neglect of the range, or tq, in 

qcot6{q) = -— + '^q^ + ... (46) 

produces 100% errors at ~ 34 MeV/c. This constraint would require 
elementary NN experiments to be performed at scattering energies < 2 MeV, 
which appears to be impractical except, potentially, in future spin rotation 
experiments in parahydrogen. 

While measuring NN scattering observables at very low energies may 
be prohibitively difficult, one might propose alternative measurements of 
near-threshold PNC observables in few-nucleon system, such as neutron spin 
rotation in ^He. Unfortunately the effective center-of-mass momentum rele- 
vant to NN partial-wave analyses is generally not determined by the external 
kinematics in such reactions, but rather by the nuclear Fermi momentum, 
typically ~ 200 MeV> m^. This is the characteristic scale of the momentum 
fiow in the exchange term, which often dominates the PNC interaction, as 
Barton's theorem restricts direct terms. There is a nice pedagogical treat- 
ment of these effects in [6] (see Sec. 6.1), where a mean-field approximation 
to V^'^'-'' is derived in a Fermi gas model, showing analytically that exchange 
terms make important contributions to every S — P amplitude, are the only 
source of sensitivity to h^., and tend to dominate couplings like /i° in light 
isoscalar nuclei. 

Consequently, we are stuck with the need for evaluating PNC observables 
in a framework that treats strong interaction distortions. This framework 
is not provided by Danilov or pionless EFT approaches in which strong in- 
teraction input is hmited to scattering lengths. Pionful EFTs could be con- 
templated, but their development for strongly interacting A > 2 systems 
has proven to be challenging due to the anomalous scales mentioned above, 
apparent from nuclear binding energies: this was the origin of Weinberg's 
proposal to distinguish irreducible graphs that involve no infrared enhance- 
ments, treating these by chiral perturbation theory, from reducible graphs, 
which involve the strong potential V in combination with infrared-enhanced 
propagators, requiring all-order summations. Such considerations have led 
to hybrid EFT strategies where wave functions are generated from potentials 
that have been tuned to the anomalous nuclear scales, but where the inter- 
actions evaluated between such wave functions are developed through EFT. 
This appears to us to be a viable strategy for modern treatments of PNC, with 
the coefficients of the five S — P operators providing a model-independent 
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description of the PNC NN interaction. Such a theoretical program would 
be timely if new data soon become available. 

This leads us to the question of what has been measured vs. what might 
be measured (or calculated) in the near future to further our understanding 
of hadronic PNC. First we consider the constraints illustrated in Fig. 3, 
where we take all calculations from [6] except as noted: 

1. The longitudinal p + p asymmetry at 13.6 [10], 15 [11], 45 [12], 221 [13] 
MeV: 

The new combined analysis described in Sec. 4 yields the S — P constraint 

A'f = A;^°-'^° + A;^"-'^» + " ^ = (4.19 ± 0.43) X 10"^ (68%c.l.) 

v6 

It also yields a weaker constraint on the — ^ -D2 amplitude which we give 
here in terms of DDH couplings 

Qph^^Xv + g^hP^Xs = -(4.4 ± 1.6) x 10"^ (68%c.l.) 

2. The longitudinal p+'^He asymmetry at 46 MeV [20] 

Al{vo^: 46 MeV) = -0.025^^^jv/il + 0.0505^/iJ + 0.017^^/iJ + 0.007^^/1° 

~ -0.00355Ai'^°"'^° - 0.00317A;'^»~'^" 
-0.00268Ao^i-'^^ - 0.00114Ai^^-'^i 
= -(3.3 ± 0.9) X 10"^ (47) 

3. The circular polarization of the 7-rays omitted in the decay of excited ^^F. 
The data from the various experiments discussed earlier can be combined to 
determine P^(^^F) < (1.2±3.9) x lO"''. As the relationship of this limit to the 
underlying weak parameters depends on a mixing ratio of known magnitude 
but unknown sign [6], one finds in terms of DDH couplings 

\P^Q^¥)\ = \-i2Qg^MNhl - llQgphl - img^hl\ < 5.1 x 10"^ 

where the nuclear matrix element has been taken from axial-charge 5-decay 
measurements. This result can be recast as an approximate constraint on 

S — P coefiicients 

|P^(^^F)| - 15.0|Ai^i-'^i + 2.42A}^°-'^°| < 5.1 x 10"^ 
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As discussed previously, there is also a significant result from lattice QCD 
that provides an important comparison for the ^^F result [48] 

/i^(LQCD connected) = (l.lO ± 0.51(stat)+J:eJ(sys)) x 10"^ 

4. The 7-ray asymmetry in ^^F can be converted to the following constraint 
on weak couplings using nuclear matrix elements calibrated by axial-charge 
P decay (see earlier discussions on associated uncertainties), 

^^(19f) = -7.00y,^jv/i^ + 12.2^,/i° + 3.65^,/iJ + 2.31yX + 2-02yX 

1.12A;^°-'^° - 0.75A;^°-'^° - 0.48Ao^i-'^i - 0.32Ai^i-'^i 

= -(7.4 ± 1.9) X 10"^ (48) 

These constrains are similar to those for p+^He, another odd-proton system. 

As we have discussed previously, the results above, which are shown 
graphically in Fig. 3, are those that we feel can be reliably interpreted. 

There is a great need for additional measurements, particularly in NN and 
few-body systems where the associated strong-interaction effects can be han- 
dled well: 

1. The photon asymmetry in np — )■ d'-y [14, 15] 

A^{np -^d + 'j) = -O.OOSOg^NNhl - 0.0005^p/iJ + 

(3.7 X 10-^)A^-'^^ 

(0.6±2.1)xl0-^ 
-1.2 ± 1.9 ±0.2) X 10"^ ^ ^ 

The SNS continuation of the program begun in [15] may be able to reach 
a sensitivity of ~ 10~^. The lattice QCD prediction of Aj ~ -1.2 x 10~^ 
and the ■'^^F results suggest that such sensitivity will be necessary to see a 
nonzero signal. 

2. The neutron spin rotation in ^He [21]. Using [67], we find 



dz 



= [-0.072^^jviv/ii - 0.115gphl + 0.039^^/1^ - 0.026gX^+ 

+ 0.026g^hl] rad/m 
~ [o.0138Ao^°-'-^° - 0.0087A;^°-'-^° + 0.0033Ao^i-'-^i 

- 0.0033A'/i-'-^^ rad/m 
= (1.7 ± 9.1 ± 1.4) X 10" W/m (50) 
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This observable is an isospin complement of Aj(^^F) and Aiipa), so it is 
unfortunate that the sensitivity necessary to obtain an important, orthogonal 
constraint has so far not been reached. 

3. The circular polarization of the photons emitted in the capture of unpo- 
larized thermal neutrons by protons [16, 17] 



Pj(np ^ d + 7) = -OMlgph^ - 0.0088^p/iJ + OmOlg^h^ 

~ -0.00012Ao^°-'-^° + 0.00105Ao^i-'^i + 0.00154A2^°-''^° 
= (1.8 ± 1.8) X 10-^ (51) 

4. The analyzing power for p+d. The existing limit was obtained at 15 MeV, 
while the theoretical estimate comes from [68]: 

Al{p + d) = -QMllg^NNK + 0.0085^p/i° + 0.0035^p/i^ 

15 MeV ^ ^ 

+0.002(7^/1° + Q.mihgX^ 
~ -O.OOlOAo^""'-^" - .0007A5^o-'-^° 

-0.0002Ao^i-'-^i - 0.0008A'/i~'-^i 
= -(0.35 ± 0.85) X 10"^ (52) 

5. The gamma-ray asymmetry 

A^{n + d^t + -f) = 0mig^NNhl-0.l2gphl + 0m6gphrho^ + 0.020gphl 

-0.027g^hl + 0.007g^hl 
- 0.0139Ai^°-'-^° - 0.0055A;^°-'-^° + 0.0037Ao^i-'-^i 
+0.0024Ai^^~'-^i - 0.0035A2^°-'-^° (53) 

A measurement of this asymmetry was reported some years ago, but is widely 
disregarded because of its size [69] . 

6. The as yet unmeasured neutron spin rotation in hydrogen 

= [-0.23g^NNhl - 0.0S2gphl - O.Oll^^/iJ - O.OQO^^/i^ 

-Omig^hl + OMlg^hl] rad/m 
= 0.015Ai^°-'-^° - O.OllAi^i-'^i + 0.016A2^"-'^»1 rad/m (54) 
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Other few-nucleon observables have been considered in the hterature (though 

we are aware of no measurements). Ideas inchidc the longitudinal asymmetry 
for "^He(n, p)^H [18]; the photon asymmetry where unpolarized neutrons are 
scattered off a polarized deuterium target, A^{nd); neutron spin rotation in 
a deuterium target [71]; and the capture of circularly polarized photons on 
deuterium 7^ — > np. 

8 Outlook and Summary 

After more than two decades during which few new results became available, 
we are beginning once again to make progress in imderstanding how the AS = 
weak interaction operates among strongly interacting nucleons. First we 
argued here that if the existing experimental results yielding nonzero values 
for PNC observables arc treated in a consistent formalism, the agreement 
among them is really quite good: The p + p, ^^F, ^*^F, and p+'^He results 
combine to suggest a ratio of isoscalar-to-isovector strengths about a factor 
of six larger than the "best value" benchmark of DDH. The tension that had 
existed among the results appears to have arisen from comparisons that did 
not utilize a common set of strong meson-nucleon vertices. 

Second, we have the first tentative confirmation of the most puzzling 
result in the field, the indication from ^^F that h]. is at least a factor of three 
below the nominal DDH "best value." The evidence supporting this result has 
come from theory, the first lattice QCD calculation of a weak meson-nucleon 
coupling. The lattice QCD value for is consistent with the ^^F upper 
bound. While the result is tentative - calculations at the physical pion mass 
with the inclusion of disconnected contributions remain to be done - large 
changes are not expected, based on current estimates of the calculation's 
statistical and systematic uncertainties. 

Third, as lattice QCD gives an near the upper bound of the ^^F band, 
the ongoing SNS experiment on n + p^d + '-f may succeed in measuring 
isovector PNC, if it reaches its precision goal of ~ 10~^. This would give us 
a direct experimental value for hi. or, equivalently, A^'^^~ 

Fourth, there is a very good prospect that could be calculated rather 
precisely from lattice QCD: this contribution has only connected pieces. This 
would determine A2'^°~ ^° and, consequently, allow one to extract from p + p 
much sharper constraints on Aq'^"^ + A^'^°~ Our calculation combining 
the 13.6, 45, and 221 MeV p + p results could be redone, without the loss of 
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sensitivity that comes from marginalizing over possible values of the unknown 

parameter Ag'^"" 

This review has also tried to emphasize the close relationship between the 
various formulations of hadronic PNC, whether based on potentials such as 
DDH [40] or pionless effective field theory, such as the calculations of Zhu et 
al. [59], Girlanda [62], and PhiUips et al. [63]. The common language is the 
set of five S-P operator coefficients that provide a model-independent param- 
eterization of the most general low-energy PNC interaction. Operationally, 
the DDH potential differs from pionless EFT interactions only through the 
inclusion of a sixth parameter that can mock up the effects of long-range 
pion-exchange in inducing — transitions; alternatively, that degree of 
freedom can be removed through a constraint, to obtain an analog of pionless 
EFT, as was done in [6]. 

With the exception of the p + p result at 221 MeV, the data we have 
utilized in this review can be represented well in calculations that employ 
the S — P operators. However, we have also stressed the importance of 
embedding that operator between realistic strong-interaction wave functions: 
momentum transfers in many of the processes of interest are characteristric 
of the Fermi momentum, and are thus well beyond the range of validity of 
Danilov or EFT treatments that hmit strong interaction input to the singlet 
and triplet scattering lengths. The importance of distorted waves for accurate 
calculations of low energy p + p, for example, has been known for many 
years. The development of pionful theories to systems with A > 2 has proven 
a great challenge, reflecting the anomalously low scale of nuclear binding 
energies and the associated difficulties with infrared enhancements. The 
kind of approach we envision succeeding in PNC studies combines a model- 
independent description of the weak interaction, in the manner of pionless 
EFT, with state-of-the-art wave functions taken from modern potentials. 
This is sometimes termed hybrid EFT. 

Finally, with lattice QCD now making contributions and with a new 
experimental program underway at the SNS, it may be time for a more co- 
ordinated theory effort on PNC. In the previous section we have described 
the dependence of PNC observables on the underlying S — P operator coeffi- 
cients. Many of the calculations used in that section, however, are dated and 
should be revisited. One purpose of this review is to set the stage for such 
a comprehensive effort. The formalism used here emphasizes the common 
features of existing PNC treatments, whether based on potentials or con- 
structed as pionless EFTs. We believe a treatment that exploits the S — P 
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operator coefficients as a "Rosetta stone" for PNC, but evaluates operators 
between wave functions that are constructed from tlie best modern poten- 
tials, applied systematically to all of the relevant few-nucleon systems,would 
address the needs of the experimental community. The simplicity of the 
S — P operator formalism would be retained, but the realism achievable with 
modern potential treatments of the strong interaction would not be sacrificed. 
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